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Abstract. We study a phenomenon occuring in various areas of quantum physics, in which an 
observable density (such as an energy density) which is classically pointwise nonnegative may 
assume arbitrarily negative expectation values after quantisation, even though the spatially in- 
tegrated density remains nonnegative. Two prominent examples which have previously been 
' studied are the energy density (in quantum field theory) and the probability flux of rightwards- 

moving particles (in quantum mechanics). However, in the quantum field context, it has been 
shown that the magnitude and space-time extension of negative energy densities are not arbi- 
trary, but restricted by relations which have come to be known as 'quantum inequalities'. In the 



X 



;_j ■ present work, we explore the extent to which such quantum inequalities hold for typical quantum 



mechanical systems. We derive quantum inequalities of two types. The first are 'kinematical' 
quantum inequalities where spatially averaged densities are shown to be bounded below. Specif- 
ically, we obtain such kinematical quantum inequalities for the current density in one spatial 
dimension (imposing constraints on the backflow phenomenon) and for the densities arising in 
Weyl-Wigner quantization. The latter quantum inequalities are direct consequences of sharp 
Carding inequalities. The second type are 'dynamical' quantum inequalities where one obtains 
bounds from below on temporally averaged densities. We derive such quantum inequalities in 
the case of the energy density in general quantum mechanical systems having suitable decay 
properties on the negative spectral axis of the total energy. 

Furthermore, we obtain explicit numerical values for the quantum inequalities on the one- 
dimensional current density, using various spatial averaging weight functions. We also improve 
the numerical value of the related 'backflow constant' previously investigated by Bracken and 
Melloy. In many cases our numerical results are controlled by rigorous error estimates. 
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1 Introduction 



The uncertainty principle lies at the root of many of the counterintuitive features of 
quantum theory. Consider, for example, a quantum mechanical particle moving in one 
dimension, whose state is a superposition of right-moving plane waves. Although the 
expectation value of (any power of) its momentum is positive, nonetheless it is possible 
for the probability flux at, say, the origin to become negative. Thus the probability of 
finding the particle in the right-hand half-line can decrease! 

We will return to this phenomenon, which has come to be known as backflow pil7ll3Uj. 
in Section 12.11 Another, related, phenomenon occurs in quantum field theory. Even 
if one starts with a classical field theory in which energy densities (as measured by all 
observers) are everywhere positive,^ one finds that the renormalised energy density of 
the quantised field can assume negative values jTTj and (in all models known to date) 
can even be made arbitrarily negative at a given spacetime point by a suitable choice 
of state. For example, the energy density between Casimir plates is computed to be 
negative; a fact indirectly supported by experiment ([HI; see the recent review jSj for an 
exhaustive list of up to date references). Various authors have suggested employing such 
effects to sustain exotic spacetime geometries containing wormholes [SI] or 'warp drive' 
bubbles |2]. Such suggestions are, however, severely constrained [221 ESI t>y the existence 
of bounds, known as quantum inequalities (QIs) or quantum weak energy inequalities 
(QWEIs) [iHllllllIllllHlliniisniiniiisllSillSl which impose limitations on the magnitude 
and duration of negative energy densities. To give an example, let {p{t))^ be the energy 
density of the free scalar field^ measured along an inertial worldline in Minkowski space. 
Then, for any real-valued smooth compactly supported the averaged energy density 



for all physically reasonable (Hadamard) states ip, where Q is a known function of poly- 
nomial growth. 

The purpose of this paper is to apply techniques developed in the field theoretic setting 
to quantum mechanical problems. In so doing we wish to draw attention to a circle of 
ideas — including sharp Carding inequalities, dynamical stability and the QWEIs — which 
eventually ought to be seen in the wider context of quantisation theory. 

We begin with a discussion of 'kinematical QIs' in Sect. |2l taking the probability 
flux as our main example. We develop bounds on spatially averaged fluxes which share 
some technical similarity with the QWEI proved by two of us for the Dirac field (see 
also 53). An important aspect of our treatment is the numerical analysis of these bounds. 
To some extent this is motivated by the recent observation of Marecki [21] that QIs may 
have observational consequences in quantum optics; it is therefore important to know 
how sharp analytically tractable bounds are. The techniques used here may also be of 
independent interest, and we give a detailed account in Sect. [31 Before that, in Sect.lSlwe 
establish a very general form of kinematical QIs arising in the Weyl-Wigner approach to 
quantizing classical systems (see, e.g., [21] as a general reference on that approach). More 

^In general relativity, one would say that the field obeyed the weak energy condition (WEC). 
^Similar statements can be made for the Maxwell, Proca and Dirac fields |32[I16[[TH] . 
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precisely, we consider the (quantized) configuration space density M" 9 x y-^ {pf{x))^ for 
normalized wave-functions ip G ^(M") which are associated with classical observables F, 
i.e. functions on phase space x R", by 



where denotes the Wigner function of ip. Even if F is everywhere nonnegative, the 
density {pf{x))^ may assume negative values owing to the indefinite sign of the Wigner 
function; in fact, we show that, under very general conditions on F, this quantity is 
unbounded above and below for arbitrary given x upon varying ip. Conversely, if F 
belongs to a certain class of symbols (in the sense of microfocal analysis [23 EH]) which 
are of second order (or lower) in the momentum variables, and if F is everywhere non- 
negative, then we establish a kinematical quantum inequality of the form 



with a suitable constant C depending on the non-negative weight-function x, but not 
on the (normalized) wave-function ijj. This is a straightforward consequence of the sharp 
Carding inequality [121 EH EE] ■ The general result will be illustrated by a direct derivation 
of a kinematical QI for the energy density. 

In Sect. |3]we focus attention on 'dynamical' QIs which bound temporal averages of 
the energy density in general quantum mechanical systems. (In fact our kinematical flux 
inequality may be regained as the special case, in which the evolution is the group of 
translations on the line.) These are conceptually much closer to the QIs which have been 
obtained in quantum field theory; in fact, the method we use to establish these dynamical 
QIs makes contact with the techniques employed in Some features of the general 

result will be illustrated by taking the harmonic oscillator as a concrete example. We 
summarise our main results in the conclusion. Sect. [HI 

2 A motivating example 
2.1 Probability backflow 

We begin with a simple example: the motion of a quantum mechanical particle in one 
dimension. Some time ago, AUcock [^ pointed out the existence of rightwards-moving 
states (in which the velocity is positive with unit probability) but for which the proba- 
bility of locating the particle in the right-hand half-line is instantaneously decreasing — a 
phenomenon known as probability back-flow. This phenomenon was subsequently studied 
in much greater detail by Bracken and Melloy (see also [SOD- To illustrate the idea, 
let us suppose the normalised state ip (as well as being square-integrable itself) has a 
continuous, square integrable first derivative. Then the corresponding probability flux at 
position X is given by 




(1.2) 




(1.3) 



Reip{x){pilj{x)) 



(2.1) 



m 
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where the momentum operator is, as usual, p = —ihd/dx and the particle has mass m. 
Now the spatial integral of the flux is 



Re {ip \ p^p) (p)^ 

j^{x) dx = = — ^ (2.2) 

mm 

and therefore yields the expected velocity. If is a. normalised right-moving wave-packet, 
it may be written by means of the Fourier transform as a superposition of right-moving 
plane waves 

/dk ^ 
^e^'^'m (2.3) 

with ipi^) = for < 0, so 

dk 



{Ph= I ^hk\m\'>0 (2.4) 
Jo 27r 

and we see that the spatially integrated flux is positive. However this does not imply that 
the flux itself is everywhere non-negative. Indeed, suppose that 

^^{k) = N'x[oM(k) [kVS-ko) , (2.5) 

where xn denotes the characteristic function of f2 and J\f = (/i;o(2 — ^/S) / {271))^^^"^ is a 
normalisation constant. One may calculate 

hkf. fl 1 \ hkl , , 

j . = — ^ ^ ~ -0.006^ , 2.6 

which is not merely negative, but can clearly be made as negative as we wish by tuning 
fco- Because the probability flux is negative at the origin, the probability of locating the 
particle in the left-half line is instantaneously increasing, thereby providing an example 
of the backflow phenomenon mentioned above. 

Backflow provides a nice illustration of the inadequacy of the phase velocity alone to 
predict the motion of a wavepacket. The three plots in Figure ^ indicate the evolution 
of the position probability density in time; although the packet moves to the right, the 
two main peaks are reshaped in such a way that net probability has passed from the 
right-hand half line to the left. The wavepacket is given by Eq. (|2.5|) at time t = with 
fco = 5, m = 1/2 and h = 1. 



2.2 A quantum inequality for the flux 

As we will see in Sect. El the backflow effect may be traced to the uncertainty principle. 
From this point of view, it is natural to seek bounds on its magnitude and extent. Bracken 
and Melloy [7j approached this question by showing that the probability P{t) of flnding 
a right-moving particle in the left-hand half-line obeys 

P(t)<P(0) + A, (2.7) 
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Figure 1: Evolution of a wavepacket, illustrating the backflow phenomenon. From left to 
right, the plots show the position probability density at times t = —0.1, t = and t = 0.1. 

for all t > 0, where the dimensionless constant A is the largest positive eigenvalue of the 
equation 



(for G L^(]R"^)). They also presented numerical evidence that A ~ 0.04. Using the nu- 
merical methods described in Sect.|31 we have recalculated this quantity to a much higher 
accuracy, although we have been unable to obtain consonant analytical error estimates. 
It turns out to be convenient to change variables to x = m^; we then consider the trun- 
cation of the resulting integral kernel to [0,X]. The maximum eigenvalue A(X) was then 
calculated for values of X ranging from 6000 to 24000, using X/2 quadrature nodes. This 
choice was based on calculations using a variety of densities for values of X around 2000 
for which X/2 nodes provide accuracy to 5 significant figures. By contrast, the largest cal- 
culation conducted in corresponds to X = 625, which reflects the increase in available 
computing power over the past decade. The resulting data may be fitted to a remarkable 
degree by the form A(X) = a + b/ \fX (as already noted by Bracken and Melloy for their 
data). Using a least squares fit to this, we obtain the estimate A = 0.03845182014 with a 
maximum percentage residual error under 4 x 10"'^%. Assuming the residual errors would 
be comparable for larger X, this suggests that A = 0.038452 to this level of precision. 
Our data points and the best-fit curve are shown in Fig. |21 

One may interpret the Bracken-Melloy bound (j2.7p as a demonstration of the transi- 
tory nature of backflow: large negative fluxes for right-moving states must be short-lived. 
Here, we present an apparently new bound, which demonstrates that such fluxes are also 
of small spatial extent, and whose proof is related to the quantum weak energy inequal- 
ities derived by two of us for the Dirac quantum field (TH] (see also [IS])- We consider 
spatially smeared quantities of the form 



which may be regarded as the instantaneous probability flux measured by a spatially 
extended detector. For any smooth, compactly supported, complex-valued function g, we 




(2.8) 




(2.9) 
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Figure 2: The least squares fit of A(X) to a + h/yX. 



will show that 



dx > 



h 



dx \g'{x)\" 



(2.10) 



for all normalised states belonging to the class M of right-moving states defined by 
^ = {ijj G L^(R) : ip{k) = for < and ip' continuous and square-integrable} . 



In fact, the conditions on both g and may be weakened slightly.^ 
Before giving the proof, let us make three observations. 



(2.11) 



1. First, we note that there is no upper hound on the smeared flux. To see this, choose 



any normalised ijj ^ ^ and let ipxix) 
moreover, 



e ip{x). We have ipx e 



so / j^^{x)f{x) dx +00 as A ^ +oo. 



Xh. 



m 



for A > 0; 
(2.12) 



2. Second, the scaling behaviour of the above bound may be investigated by replacing 
g by gx{x) = \~^^'^g{x/ X), whereupon the right-hand side of inequality (j2.10|) scales 



■^In particular, continuity of tp' may be weakened to & AC{R) fl L^(K) with ip' G i^(M) at the 
expense of augmenting some statements with the quahfication 'almost everywhere'; by an approximation 



argument it is easy to see that (|2.1(J|I holds for all g belonging to the Sobolev space W^'' 



m 
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by a factor of A^^. The limit A ^ corresponds to the unboundedness below of the 
probability flux at a point, while the limit A — oo is consistent with the fact that 
> for T/^ G ^ (because the bound vanishes more rapidly than A~^). Roughly 
speaking, our bound asserts that the magnitude of negative flux times the square of 
its spatial extent satisfies a state-independent upper bound on Thus the extent 
of backflow is limited both in space and in time. 

Note also that the bound ()2.10j) vanishes in both the classical limit /i — > and the 
limit of large mass. This differs from Bracken and Melloy's inequality ()2.7|1 in which 
the dimensionless constant A is independent of h and m. 

We remark that — again in contrast to [7^ — our result is kinematical rather than 
dynamical: no specific Hamiltonian is invoked. Here, 'kinematic' refers to the kine- 
matics of the Schrodinger representation, i.e., the (unique) regular representation of 
the Heisenberg commutation relations. 

3. Finally, on integration by parts, Eq. ()2.10|) can be reformulated as the assertion that 
for each normalised ip & the Schrodinger operator 

is positive on the space of smooth compactly supported functions g, in the sense 
that 



g{x){H^,g){x) > (2.14) 

for all such g. Although the physical significance of this reformulation is not clear, 
it can provide useful necessary conditions for a given function j{x) to be the flux of 
a right-moving state. Only if the corresponding Schrodinger operator has no bound 
states can this be the case. This can be sharpened slightly: as an illustration, 
suppose j^(,{x) is the flux of a state in ^ with j^{x) < —M on some open interval / 
of length a. Then positivity of on C^{I) implies that the Friedrichs extension 
Hm of the operator 

-AirhM on Co°°(/) C L'(/) (2.15) 



2m dx^ 



is also positive. Since the Friedrichs extension of this operator corresponds to the 
imposition of Dirichlet boundary conditions at the boundary of /, Hm has spectrum 



= ^ - iTihM (n = 1,2,3,...) (2.16) 



2ma^ 

so we deduce (from Ei > 0) that M < hir / {Sma'^) . This provides a more quantita- 
tive version of the connection between the magnitude and spatial extent of negative 
fluxes. Similar ideas have been employed in the context of quantum weak energy 
inequalities jlH] to cast light on the 'quantum interest conjecture' of Ford and Ro- 
man 
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We now establish the quantum inequahty (|2.1(Jj) . It is sufficient to prove this for the 
case in which g G is real- valued. Setting f{x) = g{x)^ and writing Mf for the 

multiplication operator {Mfip){x) = f{x)4'{x), we have 



j^{x)f{x)dx = —ReiijjlMfpip) 

= ^Re ((^ I MgpMg^) + I Mg[Mg,p]ij)) 

m 



h f dk 

m 27r 



k 



Mgij{k) 



(2.17) 



where we have used the fact that Re{i/j \ Mg[Mg,p]ip) = Reih{4' \ Mggiip) = 0. We 
therefore may obtain a bound by estimating the portion of this integral arising from 
k<0: 



h f'^ dk 
ji,{x)f{x)dx > — I —k 



Mgi^ik) 



h 



dk 



k 



Mg^{-k) 



(2.18) 



By the convolution theorem 

M^{k) 



dk' 
2^ 



^ik')gik - k') , 



(2.19) 



where the restriction to k' G is permissible for ip & M. Now a straightforward 
application of Cauchy-Schwarz gives 



Mgtlji-k) 



< 



dk' 
2^ 



'm2 



g{k + k')\ 



where we have also used |^(— A;)p = |^(/;;)P (since g is real) and 
in (j2.18|) . we now calculate 



(2.20) 
1. Substituting 



[Mx)fix)dx > -- r^k\g{k+k')\' 

J mJo 2tt Jo 2tt 



'0 ^" JO 

h r du 



m 

h_ 

m 

h_ 

m 







i2ny 



\9W 



dk k 



du 



ire 



y\9in)\' 



du 



2^ I^HI^ 



57rm 



167r 



dx \g'{x) 



(2.21) 



where we have changed variables from {k, k') to {u, k) with u = k + k', used evenness of 
\'g{u) \ and Parseval's theorem. This completes the proof of the quantum inequality ()2.10|) . 
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The later stages of this argument may be rephrased as follows. The inequality (|2.18|) 
asserts that 

[ j^{x)f{x)dx>--\\T^' (2.22) 
J m 

where the operator T acts on L^(R+, (iA;/(27r)) by 

/dk' 
— Vkgi~k - k'Mk') (2.23) 

and is easily seen to be Hilbert-Schmidt. Varying over normalised ip, the right-hand side 
of Eq. (j2.22p is bounded below by — ||T||^, where ||T|| denotes the operator norm of T. 
This leads to the bounds 

/ Mx)fix) dx > --lirf > --llTllli.s. , (2.24) 
J m m 

where the last inequality holds because the Hilbert-Schmidt norm ||T||h.s. dominates the 
operator norm. The calculation in ()2.2H1 in fact precisely computes this final bound. 

To summarise, we have seen that, even for a right-moving state ip & the flux 
need not be pointwise nonnegative; moreover by tuning the state, one may arrange the 
normalised flux j^{x) to be as negative as one likes at a given fixed x. However, weighted 
spatial averages of the flux are bounded below in terms of the weight function alone. This 
condition may be reformulated as asserting the positivity of the Hamiltonian for a particle 
moving in a potential given (up to constants) by the probability flux of any state in 

2.3 Numerical results and sharper bounds 

We illustrate our bound by reference to four weight functions: a Gaussian, a squared 
Lorentzian and two compactly supported weights which we call the truncated cosine and 
the smoothed truncated cosine. (Neither of the compactly supported weights are C°°, but 
they have sufficient smoothness for the above argument to hold; see footnote El above.) 
Our weight functions are summarised in Table ^ along with the corresponding bound 
arising from Eq. (j2.1(jp . In each case, f\ has unit integral, the parameter A controls the 
sampling width and gx{x) = a/ fx (x) . For later reference we have also given the Fourier 
transforms of fx and gx- 

We wish to compare the above bound with two sharper (but less analytically tractable) 
bounds: the bound arising from the first inequality in ()2.24j) and a direct numerical 
estimate of the infimum of the integrated fiux. In the first case, we are required to find 
the operator norm of T. Our numerical approach proceeds by first truncating the kernel 
to an interval [0, K] — we are able to estimate the error incurred here by using bounds 
obtained from the Hilbert-Schmidt norm — and applying a numerical quadrature scheme 
due originally to Fredholm (see e.g.. Sec. 4.1 of jlj or Chapter 4 of [10]) to the truncated 
kernel. This leads to a matrix whose eigenvalues approximate those of the truncated kernel 
and hence the original operator, and which can be computed using standard numerical 
packages. Full details, including a discussion of error estimates, are given in Sect. |31 This 
leads to quantum inequalities 

/hC 
3i,i^)fxix)dx> -, (2.25) 
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vjrctLio&iclii 


Squared Lorentzian 




OiliiJiJ tiled Li Llil*^CltCLl *^iJk5lilC 


fx 


(A0F)-ie-(^/^)' 


2A37r-V(x2 + A2)2 


^?(A- |x|)/Acos(x7r/(2A)) 


4^?(A- |x|)/(3A)cos(x7r/(2A))^ 


h 


g-(A«)V4 


(l + A|A;|)e-^l^l 


sm(A/c) 


47r sm(A/c) 




AA;(7r2 - A;2A2) 

V / 


AA;(A;2A2-47r2)(A;2A2-7r2) 




^^l/4g-(«A)V2 


V2A7re-^l*^l 


^i^\f\ cos(AA;) 


27r^ sin(AA;) 


7r2 - 4A;2A2 


A;V3A(7r2 - A;2A2) 


QI bound 


-V(167rmA2) 


-V(167rmA2) 


-/i7r/(32mA2) 


-/i7r/(24mA2) 




-0.01989436788 


-0.01989436788 


-0.09817477044 


-0.1308996939 



Table 1: Compendium of sampling functions considered. 



Gaussian 
K i2{K) 


Squared Lorentzian 

K ii{K) 


Trunc. cosine 
K i2{K) 


Smoothed trunc. cosine 
K n{K) 


10 -0.0048295212087 
20 -0.0048295668511 
30 -0.0048295668517 
40 

50 
60 


30 -0.002980544308 

40 

50 

60 

70 

80 


2000 -0.029012801924 
2200 -0.029012804495 
2400 -0.029012806174 

2600 -0.029012807318 
2800 -0.029012808114 
3000 -0.029012808686 


140 -0.036095566956 
160 -0.036095567038 
180 -0.036095567056 
200 -0.036095567060 
220 -0.036095567061 
360 



Table 2: Numerical estimation of infer (J) for various kernels. 



where the constant C depends on the particular weight function used and is given: 



Kernel 


C 


Accuracy 


Improvement 


Gaussian 


0.01958128485 


10 S.F. 


1.6% 


Squared Lorentzian 


(167r)-i 


Exact 





Trunc. cosine 


0.08463957004 


2 S.F.? 


16% 


Smooth trunc. cosine 


0.125047838 


At least 3 S.F. 


4.5% 



Note that we were only able to obtain fairly weak error bounds in the truncated cosine 
case. In each case, the improvement on the analytical bound is relatively small. We 
may interpret these results as showing that T = R + S where R has rank 1 and the 
Hilbert-Schmidt norm of S is small relative to that of T. This is most apparent in 
the squared Lorentzian case, in which T is itself exactly rank 1 and no improvement is 
obtained by using the operator norm. It would be interesting to understand the origin of 
this apparently general phenomenon. 

Our second numerical calculation aims to compute the infimum of the spectrum of the 
unbounded integral operator 

Pj r°° rik' ^ 
{Jvm = - —{k + k')h{k - k'Mk') . (2.26) 
2 Jq 2tt 

We proceed by truncating the kernel to the interval [0,ii"], computing the minimum 
eigenvalue ^{K) using sufficiently many quadrature points to obtain machine precision. 
We then increase K until convergence of n{K) is obtained, again to machine precision. 
Our results are given in Table |2l in which we give /i(-ft') in units of h/{m\^). Blank 
entries indicate that the computed value was identical to the last printed number in that 
column. The density of quadrature points used (per unit K) was 5 for the Gaussian, 
1 for the truncated cosine, and 5 for the smoothed truncated cosine, although higher 
densities were also used as a numerical check (40, 2, and 10 respectively). The results for 
the squared Lorentzian were rather slower to converge as the density increased (perhaps 
because the kernel fails to be everywhere smooth) and were computed using a density of 
60. For K < 80, a density of 70 was used as a check. 

To summarise, we have seen that a) the limitations of our flux QI do not lie in the 
estimation of an operator norm by a Hilbert-Schmidt norm, but rather in the earlier stages 
of the derivation (probably the estimate ()2.18j) ): b) the overall scope for improvement on 
our flux QI is roughly a factor of between 3 and 7 (in our examples), and it is clear that 
the sharp bound is not simply a multiple of our bound ()2.10p (in contrast to the situation 
for two-dimensional massless quantum fields [ZHE])- 

3 The Wigner function and Kinematical Quantum 
Inequalities 

It is worth emphasising that phenomena similar to those presented above arise naturally in 
the context of Weyl quantisation, in which the phase space aspect of quantum mechanics 
is brought to the fore. In our discussion we will consider the phase space to be R" x 
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(see, e.g., j2H] for Weyl quantisation on manifolds). We recall that the central object in 
this approach is the Wigner function defined on phase space by 

W4x,p) = (^0 J <rye^'^y'^i,{x + y)i^{x - y) , (3.1) 

where ip ^ -L^(M") is the corresponding normalised quantum mechanical state vector. The 
classical analogue of W^{x, p) would be a probability distribution on phase space; as is well 
known, however, is not itself a probability distribution because it is not guaranteed 
to be everywhere nonnegative. This has important consequences for observables obtained 
via Weyl quantisation, which proceeds as follows. 

Given an observable on the classical phase space, i.e., a smooth function^ F : M" x 
M" — > M, Weyl quantisation defines an operator F^, whose expectation values are given 
(for normalised ip) by 

{F^h = I ^^^F{x,p)W4x,p) . (3.2) 

The action of this operator may be written in the form 

-j^^nix + y]/2,p)e^(^-^)-^/V(y) . (3.3) 

Let us note that this procedure also yields a natural definition for the quantum me- 
chanical density associated with a classical observable. Namely, setting 

{pFix))^ = J ^F{x,p)W^{x,p) , (3.4) 

it is clear that the spatial integral of {pf{x))^ yields the expectation value (F^)^ for all 
F and ip (modulo domain questions^). 

Now, because the Wigner function need not be everywhere positive, we see that the 
Weyl quantisation of a nonnegative classical observable may assume negative expectation 
values. This situation is exacerbated for the densities defined above (see statement (II) 
below). However, we will show that kinematical quantum inequalities may be derived, 
under certain conditions. Indeed, these bounds are obtained as applications of the so- 
called sharp Garding inequalities in the theory of pseudodifferential operators jSElElEZI- 
It is interesting to note that Fefferman and Phong, to whom the most general sharp 
Garding results are due, were guided by intuition arising from quantum mechanics: in 
particular, the uncertainty principle. 

We begin by specifying more precisely the class of classical observables. For m G N, the 
symbol class S"^ (often denoted more precisely as S'™q) is defined to be the set of smooth 
functions F : M"' x MJ^ C such that, for each compact C M" and n-dimensional 
multi-indices a, (3, there exists a constant CK,a,f3 such that 

mDiF) {x,p)\ < Ck,.A^ + bl)™-'^' (3.5) 

^Precise growth conditions will be specified below. 

^For example, this will certainly hold for F of polynomially bounded growth and ■0 belonging to the 
Schwartz class. 
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for all {x,p) G K X R" (see, e.g., j2ni EHI for multi-index notation). By 5"^^, we denote 
the set of F G S"^ admitting a (unique) decomposition F = Fp^ + Fsub such that the 
principal symbol Fp^ belongs to S"^, is homogeneous of degree m in momentum, i.e., 
Fpr(x, Xp) = X'^Fpf{x,p) for all {x,p) G and A G M, and is nonzero except at vanishing 
momentum, while the sub-principal symbol F^ub belongs to S"^~^. 

For F G S"", the Weyl quantisation F^ is a continuous linear map from C^(M") to 
C°^(M"), so Eq. (jSI21) holds for all normalised ip G C^(M"). The density {ppix))^ is in 
fact defined (and indeed smooth in x) for all ip belonging to the Schwartz class ^(M"); 
however, it is only guaranteed to be integrable for ip G C^(M"'). 

Below, we will establish the following observations: 

(I) Suppose F G S'l^^ is real, for some m > 1. Then, for each x, the density {pf{x))^ 
is unbounded from both above and below as ip varies in C^(R"') with HV^IlLa = 1. 

(II) Suppose F G S"^ is nonnegative, F{x,p) > for all {x,p) and let x ^ C^(]R") be 
nonnegative. Then there exists a constant C > 0, depending on F and such that 

j (rxx{x){pF{x))^ > -C (3.6) 
for all Ip G with \ \iP\\l2 = 1. 

To establish (I) we may assume without loss of generality that x = 0, for which we 
have 

for normalised ip G C^(M"). Setting ipx{x) = \~"'^'^ip{x/ X), (A > 0) and making the 
obvious change of variables, 

Xh) y (3.8) 

so, bearing in mind that 

|F(0,p/A) - Fp,{0,p/X)\ < C{1 + \p/X\r-' (3.9) 
by definition of S'^^ and Eq. ()3.5|) . we obtain 

h) J 7|^^P^(0'P)^'^''^'^^(-^) (3-10) 

as A — s> 0+. It now remains to show that the right-hand side of this expression attains 
values of both signs as ip varies in C^(M"). To this end, assume (without loss) that 
-^pr(0,p) depends nontrivially on the first coordinate, pi, of p. Integrating by parts, the 
right-hand side of ()3.1U|) in the form Py{ip{y)ip{—y))\y=Q, where Py is a homogeneous linear 
partial differential operator (in y) of order m with (possibly complex) constant coefficients 
Cq. We now consider ip of the form 
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where x ^ C^(]R") is equal to unity in a neighbourhood of the origin and / G C^(R). 
For such if) we have 

Py{W)i^{-y))\y=o = QyA7bji)fi-yi))\y.=o (3.12) 

for some ordinary differential operator 

Q.. = Ec.H)'^ (3.13) 
k=o ^1 

of order 1 < q < m with constant real coefficients. (That Qy^ is of order at least one is 
a consequence of our assumption that Fpr(0,p) depends nontrivially on pi, reality of the 
Cfc holds because the right-hand side of Eq. ()3.10p is manifestly real for all ip e C^(M").) 
We now choose / so that /(O) = 1 and f^''\0) = for 1 < A; < g — 1. Then by Leibniz' 
formula, 

PyiWm~y))\y=o = c,i-ty (7^(0) + (-i)v^^no)) +co. (3.i4) 

It is now obvious that / may be chosen so that the right-hand side of this expression 
adopts values of both signs, completing the argument. 

Statement (11) is straightforward: because x ^ (^^(M") and F G 5*^, the symbol x^ 
obeys uniform bounds 

\{D:dP^xF){x,p)\ < Q^(l + bD-^-l^l (3.15) 

for all {x,p) G M" X M", so the sharp Garding inequality (Corollary 18.6.11 in |2Z] with 
5 = 0, p = 1; see also Eq. (18.1.1)" therein) entails the existence of a constant C such 
that 

J d^xxix){ppix))^ = {{xFU^ > -C (3.16) 

for all normalised G ^(M"). This is the required kinematical quantum inequality. 
We now give two examples to illustrate the above ideas. 

Example 1: Consider a classical Hamiltonian 

2 

H{x,p) = ^ + V{x) (3.17) 
on X M'^ with V G C°°(]R"). The Hamiltonian density obtained from Eq. (jSSl) is 

ipHix)h = ^ (|V^(x)P - ReV<^(AV^)(x)) + V{xmx)\' , (3.18) 

where A = is the Laplacian. Clearly {ph{x))^ may be made arbitrarily negative as 
ip varies in C^(R") by arranging that Vip{x) = 0, tlj{x)Ailj{x) > and then — as in (1) 
above — scaling ijj about x, introducing 

^A(2/) = A-"/2^((y-x)/A) (3.19) 
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for which 



{Ph{x))^^ = -A-("+2)_ReV;(x)(AV^)(x) + \--V{xmx)\^ . (3.20) 

As in the proof of (I), the subprincipal symbol drops out in the hmit A — > 0+, so 

{Ph{x)).^^ -oo 

Since H E S"^ we aheady know that a kinematical quantum inequahty exists. However 
it is instructive to give a direct argument for this, which also yields an explicit bound. To 
this end, we note that, for any nonnegative x ^ C^(M") and normalised ip G o$^(]R"), 

[ rxxix){pHix))^ = -^{p.tP I M^p.tP) + ^{tp I (M^ +p2Mj^) + {tp I M^vi^) 

J '±711 OTJI 

(3.21) 

where pi = —ihVi and {Mfilj){x) = f{x)ip{x) is the operator of multiplication by /. Now 
[M^,p]=thM^^^, so 

M^p2 + p^M^ = 2piM^pi + ih[M^7^^,pi] = 2piM^pi - h^M^x (3.22) 

and hence 

j rxx{x){pH{x))^ = ^{p^^ I M^p.tP) + I Ml^) (3.23) 

where ^ 

L{x) = -^{^X){x) + V{x)x{x) . (3.24) 

Since the first term in ()3.23|) is nonnegative, we obtain the quantum inequality 

J rxxix){pH{x)).^ > ig^ (^-^iAx){x) + Vix)xix)^ , (3.25) 

for all nonnegative x ^ C^(M") and G ^(M"). We note as a curiosity the appearance 
of a Schrodinger operator applied to the weight x (rather than the state ip). In the case 
of a nonnegative potential, we may obtain a QI (slightly weaker than that given above) 
in the form 

/ rxx{x){pH{x))^ > ] inf {Hx){x) . (3.26) 

Example 2: We now show that (II) allows us to deduce the existence of a kinematical flux 
QI on rightwards moving states. Let / be a nonnegative smooth compactly supported 
function. The averaged probability flux j^(/) is easily seen to be the expectation in 
state ip of the Weyl quantisation j(/) of f{x)p/m, which is a (first order) element of the 
symbol class S"^ (but is of course negative for p < 0). Now let ri{p) be smooth, vanishing 
for p < and equal to p for p greater than some po, and set F{x,p) = f{x)rj{j))/m. 
Then the quantisation differs from j(/) on ^ fl S^{EP) only by a bounded operator. 
Accordingly (II) entails that j^^if) is bounded below for normalised ip E ^ (1 ^(M). Of 
course, this argument does not determine the magnitude of the bound, in contrast to the 
direct approach of Sect. 12.21 
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We should like to remark that in jB] there appears a result which is complementary to 
ours; in that reference, the authors consider the one-dimensional case n = 1 and show that 
there is a V'-independent bound below (and above) on the integral of the Wigner function 
over elliptic sub-regions of the phase-plane which is much sharper than that implied by 
the a priori uniform bounds on the Wigner function. This is again an effect of averaging, 
this time over a region of finite extension in both x- and p-space. It would be interesting 
to see if this result can be generalized to higher dimensions through a generalization of 
(II) to a more general class of symbols; however, it is not at all clear that this can be 
accomplished as it apparently goes beyond the scope of sharp Carding inequalities. 

4 Dynamical quantum inequalities 

In this section we turn to a different type of QI, which is closer to those studied in quantum 
field theory. The focus here is on time-averages of the energy density at a fixed spatial 
point: we will refer to the QI bounds obtained as dynamical quantum inequalities. To 
keep the discussion fairly general, we assume that configuration space M is a topological 
space carrying a measure z/, so that the state space is = L^(M, dz/). The dynamics 
is assumed to be generated by a self-adjoint Hamiltonian H which is defined on a dense 
domain in J^. Each normalised state ip belonging to the domain of H then determines 
both a position probability density {p{t,x))^ and a Hamiltonian density {h{t,x))^ by 

{p{t,x))^ = IM-^ (4.1) 
{h{t,x))^ = ReMx){HtPt){x) (4.2) 

where tpt = e'^^^^'^tp. This definition of the energy density differs from that employed 
in Sect. El note that we are not assuming in this section that H is the quantisation of a 
classical observable, in which case the above would appear to be the most natural defini- 
tion. In particular, both the quantities defined are integrable with respect to the measure 
(iz/(x) for each t G M, with integrals equal to unity and {H)^ respectively. However, we 
will be interested mainly in time averages of these quantities at some fixed point x G M. 
In so doing, we immediately encounter the problem that it does not generally make sense 
to speak of the value of an L^-function at a point. ^ To avoid this, we introduce the spaces 
Ji^k = -D((l + W^Y^"^) and assume that, for some A; > 0, each element in Mk should be 
(almost everywhere equal to) a continuous function and that for each x G M there is a 
vector rj-j. G J'^ such that 

ip{x) = {7],\pk{H)iP) V^G^fc. (4.3) 

Here, il){x) means the value at x of the continuous function to which ip is almost ev- 
erywhere equal, and we have written Pk{E) = (1 + E'^Y^'^. Therefore, the functional 
^ {Vxl Pk{H)ip) on coincides with the 5-distribution concentrated at x, so that 
formally [as it is not an element of Jif] pk{H)*r]x is the 5-distribution. In practice, these 

^Elements of the space L^{M,dij) are really equivalence classes of functions agreeing almost every- 
where. 
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assumptions are fairly mild: in particular, for the case in which H is minus the Lapla- 
cian on some manifold they are simply a transcription of the content of Sobolev's lemma. 
We remark that the necessary regularity in quantum field theoretic quantum inequalities 
is obtained by restricting to the class of Hadamard states, which would correspond to 
Jifoo = OkeN'^k in the present context. 

It now makes sense to define the position and Hamiltonian densities as 

= \ piH)^Pt)\' (4.4) 
{h{t,x))^ = Re{p{H)^t\Vx){v.\Pk{H)H^t) (4.5) 

for normalised states G ^k+i- Furthermore, one may easily check (using Cauchy- 
Schwarz) that these quantities are bounded in t, so the time-averaged quantities {px{,f))^ 
and (Kif))^ given by 

{pAf))i^= [ dtf{t){p{t,x))^ (4.6) 



and the analogous equation for (h^^f))^ are well-defined for any smooth compactly sup- 
ported function /. 

From now on, we denote the spectral measure of H by dPg. (In the case where H 
may be diagonalised by a basis of orthogonal eigenvectors 0„ with simple eigenvalues En, 

I I PEV)fiE) = 5^ I 0„)(0„ I v)fiEn) 

= I dEY,KE - En){iJ I I '^)f{E)- (4.7) 

more generally, the projection- valued measure allows for the case of varying-even infinite- 
multiplicities and for both continuous and discrete spectrum.) While there is some am- 
biguity in choosing k and rjx such that Eq. ()4.H|1 holds, the measure on M defined by 

/i,(A)= / {r^x\dPE'nx)Pk{Ef (4.8) 

for bounded Borel sets A has an independent meaning. In fact, Px{^) is simply the 
diagonal P/\{x,x) of the integral kernel of the spectral projection P^ of H on A,^ given 
by 

Px{A)= Yl (4.9) 

n:E„£A 

if H has purely discrete spectrum. Below, it will occasionally be useful to consider the 
corresponding measure arising from self-adjoint operators other than H; in these cases, 
we will write px to denote the operator H' involved. Finally, since < /ix(A) < 
WrjxW^ sup^^/^Pk{E)'^, we see that /ix is polynomially bounded. 

After these preliminaries, we come to the statement of our dynamical quantum in- 
equalities. 

''Since, for any V' G we have PaV" ^ D{p{H)), it follows that {PAi^){x) = J di'{y) PA{x,y)^{y) 
where PA{x,y) = {pk{H)PA'nx){y) is continuous in y. This last quantity may easily be expressed as 
{p{H)PaVx I Pk{H)PAVy), so in particular, Pa{x,x) = \\pk{H)PAVx\\^ = Ma;(A). 
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(III) Let g be any real- valued, compactly supported function on M and set f = g^- Then 
given real numbers a < b, the inequalities 

b{pAm + J ^Q+{u)\g{u)\'>{KU)h>a{p,U)h-j ^Q-(«)l?Wr 

(4.10) 

hold for all normalised if) G P\^a,b]'^i where 

Q-{u) = I dpx{E){hu + a — E}j^ 

J[a,b] 

Q+iu) = [ dp^{E){hu~h + E}+ (4.11) 

J[a,h] 

are nonnegative, monotone increasing and polynomially bounded in m, and we 
have used the notation {A}+ = max{0, A}. (Similarly, we will write {A}_ = 
min{0. A}.) Moreover, the first (resp., second) inequality in ()4.1()|1 also holds for 
all ip e P(_oo,fe]^^+i (resp., P[a,oo)'^k+i) provided the integration range in ()4.1H) is 
replaced by (—00,6] (resp., [a, 00)). 

(IV) Suppose dpx{E){l + \E\) < oc and let g be as in (III). Then, for any fixed c G M, 
the inequality 

{hAf))i.>c{pAf))^- I ^S{H-cl;u)\g{u)\' (4.12) 

holds for all ip G J^k+i, where 

S{H; u) = J dpi^\E) {hu - E}+ (4.13) 

is nonnegative, monotone increasing and polynomially bounded in u. (There is of 
course a dual statement, for the case dpx{E){l + \E\) < 00.) 

Before proceeding to the proof of these statements, we illustrate them by drawing 
some consequences. The interpretation of (III) is that a state with energy between a and 
b has an averaged energy density between a and b, suitably weighted by the averaged 
position probability density, modulo a certain latitude bounded by quantum inequalities. 
Replacing g by g\{t) = \~^/'^g{t/ A), we may consider the two regimes A — > 0"*", representing 
tightly peaked averages, and X ^ 00, which represents widely spread averages. In the 
former case, we have 

du 2 ^/^x([a', &]) du ,,2 



^^QAu)\9xiu)\'- ' " / ^H?Mr (4.14) 

provided /ia;([a,6]) > (failing which the left-hand side vanishes identically). Thus the 
latitude afforded by the quantum inequality bound grows as the sampling becomes more 
tightly peaked. As A — > 00, the QI latitude tends to zero and one may show that 

b > hmsup ^M4^ > liminf ^^4^ > a (4.15) 
A^oo {px{gi)}-4> {Px{9i))i, 
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for all ^ P[a,b]'^y provided {p{t,x))^ is nonzero for some t. This ergodic result shows 
that the spatial and temporal averages of energy densities obey related constraints. 

As a second illustration, consider (IV) in the case where H has a discrete spectrum 
{En} with corresponding orthonormal eigenfunctions {(pn}, and satisfying the integrability 
condition on fi^ (for example, H might be semibounded). Then, in the case c = 0, 

^ ^ E„<u 

= -5^«n|0„(x)p (4.16) 



a„= / du\g{u)Whu- En) . (4.17) 



where 



These formulae may be used to compare the relative ease of obtaining negative energy 
densities at different spatial locations. For example, the eigenfunctions 0„ of the harmonic 
oscillator H = /{2m) + Imtu^x^ on L^(]R) obey the following bounds (cf. the Appendix 
to Sec. V.3 in [33]): For any j G No there exists Cj > and rj G No such that 

sup 1(1 + x^)Mx)\ < + n)"-^ (4.18) 

for all n G Nq. Thus, for all normalised ijj G =5^(M) (in fact, for all ijj in a. considerably 
larger domain) we have 

oo 

{K{9')U > E ""(1 + ^)^^ • (4.19) 

' ' n.=U 

In this case, it is clear that — for a fixed sampling function g — the form a rapidly 
decaying sequence and so the sum converges for any j. Thus we have shown that the 
state-independent bound on energy density is itself a rapidly decaying function of x. 
It is therefore generally easier to maintain negative energy densities near the classical 
equilibrium point rather than far away. 

Finally, consider (III) for the case H = —id/dx on Jif = L^(]R) and a particle of 
mass m. In this instance, the spaces Jifk coincide with the Sobolev spaces W^''^(M) and 
Sobolev's lemma permits us to take k > 1/2. Then the dynamical evolution amounts 
to spatial translation and the averaged Hamiltonian density is related to the spatially 
averaged probability flux by 

{K{f))^ = mMh) (4.20) 

where f^it) = f{x — t). Moreover, the measure fi^^ is easily seen to be given by /^^^(A) = 
\A\/{2Trh), where | • | denotes the usual Lebesgue measure. Then the second inequal- 
ity in (III) may easily be checked to reproduce the flux inequality ()2.10|) for all ip G 

P[o,oo)W''+''\R). 

The remainder of this section is devoted to the proof of statements (III) and (IV). These 
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assertions are based upon two facts, which will be proved below. First, for any c G M and 
normalised ip G ^k+i-, one may show that 

{KU))^-c{pM))^ = I ^{e-c)\{^\pk{H)g{h-'[el-H])7^,)\\ (4.21) 

Second, if A is a Borel set then 

\ Pk{H)g{h-'[et - H])r].,)\' < [ dp,{E) \g{[e - E]/h)\\ (4.22) 

J A 

for all normalised ip G P^Jif. Putting these together, we obtain 

{h.im-cip^im > J ^{e-c}_Jjp^{E)\g{[e-E]/h)f 



du 

— |^(m)|' / dp,{E){hu + E-c}. 

du 
2^ 



A 



giu)\' / dp^{E){hu + c-E}+, (4.23) 



A 



for all normalised ip G Pa^, where we have made the change of variables u —u 
and exploited the fact that |^(m)P is even (because g is real- valued) . The interchange of 
variables employed in the first step is justified provided the inner integral in the last line 
of ()4.23p is polynomially bounded in u. 

To obtain the second inequality in (III), we set c = a and A = [a, h] and observe that 
the inner integral in Eq. ()4.23p reduces to Q-{u) and is polynomially bounded because 
Px is. The inequality clearly remains true for ip G P[a,oo)'^k+i with A = [a, oo). To obtain 
(IV), we set A = M and observe that the integrability condition dpx{E){l + \E\) < oo 
and polynomial boundedness of fi^ guarantee that S{H — cl, u) exists and is polynomially 
bounded. Inequality ()4.12j) follows from the above on observing that j dpi^\E) F{E — 

c) = Jdpi''-^^\E)F{E). 

To obtain the first inequality in (III) and the dual statement to (IV), one argues in 
an analogous fashion from the calculation 



{K{f))i.-c{pM))^ < I ^{e-c}^ I dp,{E)\g{[e-E]/h)\ 



2 



A 



'^'^ \giu)\' I dp,{E) {hu + E- c}+ , (4.24) 



27r 

which holds for all normalised ip G P^^- 

It remains to prove the two facts presented as Eqs. (j4.2H) and (j4.22j) above. First, 
observe that for any normalised ^ '^k, {Px{f))^ may be expressed as 

(pM))^ = I dtfit) [ d{ij I PeV.) [ d{r], \ Pij,^)e*(^-^')*/V(^)Pfc(^') (4-25) 



by the functional calculus. Performing the t integral first (which is legitimate since / is 
smooth and compactly supported) we obtain 

(Px(/))^ = jd{^\ PeV.) J d{v, I PE'i^)mE' - E]/h)pk{E)pk{E') . (4.26) 

20 



Since f = g^, the convolution theorem may be used to write 

f{[E' - E]/h) = J ^g{[E' - eynWlE^W) (4-27) 



using the fact that ^(A) = ^(— A) since g is real- valued. Substituting in ()4.2(ij] . and again 
rearranging the order of integration, we obtain 

de 
de 



{pAf))^ = i Pe^I.) / I PE'^mE' - e]/n)g{[E - e]/h)pk{E)p,{E') 



d{^ I PEr].)pk{E)g{\E-t]/h) 

= I PkiH)g{h-'[et - H])r^,)\' (4.28) 

To treat {h^^f))^ for normahsed ip G J^k+i, we write 

{K{f)h = \j dtfit) J d{i, I PeV,) j d{r], I PE'^P)e'^''-^'^'^\E + E')pkiE)pkiE') , 

(4.29) 

by functional calculus and use the identity 

^^±^/([E' - E]/h) = I ^eg{[E' - e]/h)g{[E - e]/h) (4.30) 

in place of the convolution theorem. (See JH] and J3] for proofs of this identity.) 
By a derivation analogous to that used for {px{f))^ we then obtain 

(/^.(/)U = j ^e\{,lj \ pk{H)g{h-\el - H])r^x)\^ (4.31) 

and Eq. ()4.21|) follows from this equation and ()4.28|) . 

The second assertion, Eq. ()4.22p . is proved by noting that 

|(^ I pk{H)g{h-\el - H])t^x)\^ < \\P^pk{H)g{h-^[et - H])r]x\\'' (4.32) 

using if) = P/^if) and the Cauchy-Schwarz inequality (with H^/'H = 1). The right-hand side 
may be written 



d{r,x I PEV.)Pk{Ey\g{\e - E]/h)\' = / d^^x{E)\g{\e - E]/h)\' (4.33) 

A JA 

which completes the derivation of Eq. (j4.22j) . 

5 Numerical Details 

In this section we provide more details on the numerical methods employed in Sect.|21and 
discuss rigorous error estimates on the numerical errors. The basic numerical method is 
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easily explained (see, e.g., Sec 4.1 of j3] or Chapter 4 of jl^). Suppose T is an integral 
operator on L^(M+,(iA;) with kernel G, i.e., 

POD 

{Ti)){k)= dk'G{k,k')i){k'). (5.1) 



To handle this numerically, we first truncate the kernel to [0, K] x [0, K] for some i^' > 0, 
which amounts to studying a compression of T. Provided that the required properties 
of T are, for sufficiently large well-approximated by the corresponding properties of 
Tk restricted to L^(0, K), we proceed to approximate this restricted operator by a matrix. 
To do this, we suppose that (^j)jLo ^^nd {wj)^^^ are the nodes and weights for a suitable 
quadrature method on [0,/^], and define the (A^ + l)-square matrix A = {Ajk)ff,^Q with 

{j,k) entry Ajk = w^J'^'uiJ'^G\^K{ij-,ik)- The relevant computations are performed on A 
and, if and K are sufficiently large, this will provide a numerical approximation to the 
required quantity. 

This technique was applied to the Bracken-Melloy kernel as described in Sect. 12.21 
In that case, we were unable to derive useful error estimates. However, the operator 
norm calculations of Sect. 12. HI are more controlled, as we now describe. The problem is 
to estimate the squared operator norms of the family of integral operators T\ (A > 0) 
defined in the above fashion® with kernel 

G^{k,k') = ^Vkg^{-k~k'). (5.2) 

Now the compressions Tx,x converge to Tx in the Hilbert-Schmidt norm, and therefore 
in operator norm, as K — > oo. That the matrix approximations have operator norms 
converging to ||rx,A|| as ^ oo is a consequence of the convergence of the quadrature 
formula to the integral for continuous functions. Thus our technique may be validly 
applied to this problem and it remains to control the errors inherent in the scheme for finite 
N and K. In general we have analytical control of the truncation errors (parametrised by 
K) but not the discretisation errors; we are, however, able to observe apparent convergence 
to machine precision in most cases. As a first step in our analysis of the truncation errors, 
we eliminate the parameter A from our considerations: Gx and the dilation of Gi by a 
factor of A differ only by a constant factor of A, so if truncation of Gi at K has a relative 
error of e (in either Hilbert-Schmidt or operator norm) then truncation of G\ at \K also 
has relative error e, and the associated quadrature matrices are identical apart from a 
factor of A. From now on, we will use the value A = 1 only, and write Tk for Tk^i-, etc. 

In order to estimate truncation errors, the following observations are useful. We wish 
to integrate |Gp over the region [0, oo) x [0, oo) \ [0, /T] x [0,i^]. By symmetry we can 
integrate (|G(fc, k')\'^ + \G{k' , k)\'^)/2 over the same region to obtain the same result. This 
has the advantage that 

|G(fc,fcOp + |G(A:^fc)p _ [k + k'm-k-k')\^ 

87r2 - 87r2 ^^'^^ 

which depends only on k + k'. We can exploit this by using the following decomposition 
of the quadrant [0, oo) x [0, oo) 

®In Sect. 12.31 the operators were defined on L^(M+, (ifc/(27r)). Here, we absorb the factor of (Stt)^^ 
into the kernel, which leaves the spectral data and operator norms unchanged. 
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2K 



K 



i?2 , 


i?3 







K 



2K 



We wish to evaluate the integral over i?i U i?2 U -R3, and by symmetry the integrals over 
Ri and R2 are equal. Changing into a {k + k',k — k') coordinate system we see that the 
required integral is 



47r2 



2K 



u{u- K)\g{-u)\^du + 



1 



K 



u^\g{—u)\^du . 



(5.4) 



We now consider the four sampling functions used in Sect. 12.31 Starting with the 
Gaussian kernel 

g{k) = y27ri/^e-^'/' , (5.5) 

the Hilbert-Schmidt norm of T can be found by substituting i^' = in ()5.4jl and evaluating 
the integral to give 

niH.s. = (5.6) 

which is of course an upper bound for the operator norm. 

For more precise results, we turn to the quadrature method described above. For this 
kernel, the integrals in ()5.4|1 can be evaluated explicitly to give a relative error in the 
Hilbert-Schmidt norm of 



|T-Ti 



XllH.S. 



H.S. 



;i + erf(2i\:) -2erf(i\:)) 



1/2 



(5.7) 



It can be numerically verified that the relative error falls below e = 0.5 x 10~^° (for 
ten-digit precision) at approximately K = 6.756 (this calculation requires about 25-digit 
precision) . 

The computations were performed in Maple 8 using c-panel repeated Clenshaw-Curtis 
quadrature (see Section 2.4.4 of jUIj) on the interval [0,6.9] and Maple's NAG-based 
SingularValues routine. Using 33, 65 or 129 samples with c = 1, 2 gives in each 
case the same results for the first largest two singular values: Ui = 0.1399331442, (T2 = 
0.0175697912 to 10 figures. Notice that the second singular value is very much smaller 
than the first, which means that the matrices, and hence T, can be well approximated 
by operators of rank 1. This is consistent with the operator norm, computed here to be 
.1399331442, being close to the Hilbert-Schmidt norm, l/(4v^) = .1410473959 .... This 
similarity finally justifies our use of truncation constants based on the relative error in the 
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Hilbert-Schmidt norm: the calculated value of the operator norm is certainly no larger 
than the true value, since it is the norm of a compression of T, so we have 



IT-T, 



K\ 



\T\ 



< 



IT - T.- 



A'llH.S. 



ITI 



H.S. 



ITI 



H.S. 



ITI 



< .4873572016 x 10 



-10 



.1410473959 



.4912379017 x 10 



-10 



(5. 



.1399331442 

which is still less than the target figure of 0.5 x 10~^°. 

The next kernel of interest is the squared Lorentzian; however, in this case T is a 
rank-1 operator so the Hilbert-Schmidt and operator norms coincide and there is no need 
for numerical investigation. This leaves the two compactly supported kernels. In the 
truncated cosine case, the same techniques as above lead to an error estimate of the order 
of 2% relative error in the Hilbert-Schmidt norm for K = 1100. As this is a rather weak 
estimate, we suppress the details; the numerical estimate of the squared operator norm 
(for K = 1100, N = 1024) is given in Sect. O 

Our last example is the smoothed truncated cosine, defined by 



sin(fc) 



(5.9) 



3 (7r2 - k^)k 

The relatively slow decay of this function makes the precision obtained in the Gaussian 
example impractical, but we can obtain results to at least four significant figures. In fact 
the numerical results appear to be much more precise than would be suggested by this 
error estimate. 

Maple is able explicitly to evaluate the integrals in ()5.4|1 to give a rather complicated 
formula involving the Si and Ci special functions, and from this to give the asymptotic 
formula 

for the relative error in the Hilbert-Schmidt norm. Using only the leading term, we can 
predict that truncation at about K = 732.3 should give a Hilbert-Schmidt norm relative 
error less than 0.5 x 10"'^; numerical investigation of the exact formula near this point 
confirms this value. Proceeding in the same way as for the Gaussian kernel but this time 
using the faster numerical engine in Matlab 6 to calculate the singular values, we obtain 
the following results (A^ + 1 samples, c panels): 



N 


c 


0"1 


0-2 


256 


1 


.3536210415 


.0733902393 


256 


2 


.3536211355 


.0733900951 


512 


1 


.3536210388 


.0733902259 


512 


2 


.3536210388 


.0733902259 


1024 


1 


.3536210388 


.0733902259 


1024 


2 


.3536210388 


.0733902259 



Once again, the fact that the second singular value is considerably smaller than the 
first can be used after the fact to justify the use of relative errors in the Hilbert-Schmidt 
norm (rather than in the operator norm) in choosing the truncation constant. 
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Although the error analysis only allows us to be confident of the first four figures, it 
seems likely that this figure for the operator norm is considerably more accurate than 
that. Doubling the truncation constant and using 2049 points, again with 1 and 2 panels, 
gives exactly the same results to ten figures as the two 1025-point methods above. 

The last set of calculations reported in Sect. 12.31 concern the unbounded operator J 
of Eq. ()2.26|1 . Here we have not succeeded in obtaining usable estimates of the errors 
introduced by truncation to [0,^]. However, it is nonetheless true that infcr(jR-) 
inf cr( J) as a consequence of the following arguments. 

Proposition Suppose k is an absolutely bounded kernel on L^(0, oo), and let w be a 
measurable function on (0, oo) (with respect to Lebesgue measure). Let 

D = {/ G L2(0,oo) : G L2(0,oo)}. (5.11) 

Suppose {w{x) + w{y))k{x,y) is a Hermitian function of x and y, and define an operator 
with domain D by 

poo 

{Tf){x)= {w{x)+w{y))k{x,y)fiy)dy (5.12) 
and assume that T is bounded below. For K > define the truncated operator 

{TKf){x)= f {w{x)+w{y))k{x,y)f{y)dy. (5.13) 
Jo 

Then 

lim inf a(T^) = m.ia{T) . (5.14) 

ii"— >oo 

Proof: Tk is a compression of T so inf cr(Tft-) > inf cr(T) for all K. For any £ > we can 
choose f eD with ||/|| = 1 and such that (/ | Tf) < inf ct(T) + e/2. Now 



{f\Tf) = J^ ^ {w{x)+w{y))k{x,y)f{y)dyjf{x)dx (5.15) 

and the integrand here is in L^((0, oo) x (0,oo)) by the lemma. It now follows from 
Lebesgue's dominated convergence theorem and Fubini's Theorem that, provided the 
above repeated integral can be interpreted as an integral on the measure space (0, oo) x 
(0,oo), 

(/ I 7"/) = ^ini^ ^ (^j^ {w{x) + w{y))k{x,y)f{y)dy^ f{x)dx. (5.16) 

If we let fxix) = f{x)x{Q K){x) then this tells us that {fx \ T^fx) (/ I Tf); we Iso 
have ll/xll ^ 11/11 = 1 asK ^ oo, so {fx \ T;,/K)/||/xf ^ (/ | Tf) as K ^ oo. In 
particular, for sufficiently large K we have 

^W^^""^ < (/ I Tf) + '-< inf a(T) +£ (5.17) 
II /A' II 2 

which implies that infcr(Tx) < inf o"(T) -\-e. In combination with the earlier inequality 
inf aiTx) > inf ciT), this establishes the result. 
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It remains to justify the treatment of the repeated integral as an integral on a product 
measure space. 

Lemma In the notation of the above theorem, for any f,gED, the repeated integral in 
the bilinear form 



{w{x) + w{y))k{x, y)fix)g{y)dy j dx (5.18) 

is absolutely convergent, and so can be interpreted as the integral of a function in L'^{{0, oo)x 
(0,oo)). 

Proof: We calculate 



I {w{x) + w{y))k{x, y)f{x)g{y)\dy j dx 
< I { j\\wix)\ + \wiy)\)\kix,y)\\fix)\\g{y)\dy^dx 

' POO \ 

\f{x)w{x)\j \k{x,y)\\g{y)\dy\ dx+ 

|/(x)| / \k{x,y)\\wiy)g{y)\dy]dx. (5.19) 



In the first term, the inner integral defines an function of x since k is an absolutely 
bounded kernel, and f{x)w{x) is an function of x since f E D. The first term is 
therefore finite by the Cauchy-Schwarz inequality. The second term is finite by similar 
reasoning: w{y)g{y) is an function of y since g E D, so the inner integral define an 
function of x, and / is by hypothesis. The final conclusion now follows from Tonelli's 
Theorem. 



6 Conclusion 

The main focus of this paper has been to draw attention to links between the failure 
of pointwise energy conditions in quantum field theory and a range of similar situations 
in quantum mechanics. In addition we have seen that there are links at the technical 
level between the QIs developed in quantum field theory and those obtained here in the 
quantum mechanical setting. In addition, we have made contact with the ideas and 
methods of Weyl-Wigner quantisation and sharp Garding inequalities. In conclusion, we 
briefly summarise the new results we have obtained along the way. 

First, we have seen that the backfiow phenomenon is limited in space (as well as in 
time [7j) as shown by our flux QI (j2.1Up . In particular, the magnitude of the negative flux 
times the square of its spatial extent is bounded above for all right-moving states in J^. 
We have also provided an improved numerical estimate of Bracken and Melloy's backfiow 
constant, and also given numerical evidence to support the conjecture that our flux QI 
is generally within an order of magnitude of the optimal bound (i.e., the infimum of the 
spectrum of J, given by Eq. ()2.2(j|l ). 
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Second, we have shown that similar phenomena occur for densities of observables ob- 
tained via Weyl quantisation. This is a consequence of the indefinite sign of the Wigner 
function, and therefore an expression of the uncertainty principle. Moreover, for observ- 
ables which are second order (or less) in momentum, we have seen that sharp Garding 
inequalities entail the existence of kinematic quantum inequalities. We have also obtained 
explicit bounds in the case of Schrodinger operators with smooth potentials. 

Finally, for general quantum mechanical systems describing dynamics on a topological 
measure space, we have shown that the time-averaged energy density obeys dynamical 
quantum inequalities (evolution being generated by the spatial integral of the energy 
density). For the 1-dimensional harmonic oscillator, we saw that the QI bound (for a 
given sampling function) is a Schwartz-class function: it becomes rapidly much harder 
to create sustained negative energy densities away from the classical equilibrium point. 
Moreover, we have seen that a bound on the spectral behaviour of the Hamiltonian on 
the negative spectral axis, expressed by the integrability condition on fi^ in (IV), already 
leads to dynamical QIs. This integrability condition can be viewed as a condition on the 
global dynamical stability of a quantum system, much in the sense of quantum systems 
in thermal equilibrium, where the spectral weight of the generator of the time-evolution 
(the Liouvillian) is exponentially suppressed on the negative half-axis (cf. Prop. 5.3.14 in 
[8j). This indicates again the link between (thermo) dynamical stability and dynamical 
QIs which was established in [20] and which originally motivated the introduction of QIs 
in [22]. 

All these findings corroborate the intimate connection between QIs and the fundamen- 
tal principles of quantum mechanics: the uncertainty principle and dynamical stability. 
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